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Abstract 

In this article, we calculate the B* —± r] c form-factors with the three-point QCD sum rules, 
then study the semileptonic decays B* —¥ rj c £i>t, the predictions can be confronted with the 
■ experimental data in the future at the LHCb. 
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1 Introduction 



The bottom-charm quarkonium states are of special interesting, the ground states B c and B* which 
lie below the BD, BD* , B*D, B*D* thresholds cannot annihilate into gluons, and decay weakly 
through b — > cW + , c — > sW + , cb — > W + at the quark level, furthermore, the B* mesons also have 
the radiative transitions B* — > B c ^. The B^r mesons have measurable lifetime, while the B*^ 
Q-l mesons would have widths less than a hundred KeV pQ. The semileptonic decays B^r — > J/ipi^vg, 

B+ —> J/ipe + D e were used to measure the B c lifetime and the hadronic decays Bf —> J/ipir^ 
were used to measure the B c mass in pp collisions at the energy ^fs — 1.96 TeV by the CDF 
and DO collaborations [2j El HJ H] . Now the average values are tb c = (0.45 ± 0.04) x 10~ 12 s and 
TO s c = (6-277 ± 0.006) GeV from the Particle Data Group [6]. The B* mesons have not been 
observed yet, but they are expected to be observed and their properties be studies in details at the 
Vy~) , large hadron collider (LHC). The LHC will be the world's most copious source of the b hadrons, 

and a complete spectrum of the b hadrons will be available through gluon fusion. In proton-proton 
collisions at ^/s = 14 TeV, the bb cross section is expected to be ~ 500/i& producing 10 12 bb pairs 
Q\ • in a standard year of running at the LHCb operational luminosity of 2 x 10 32 cm~ 2 sec -1 [7]. 

The semileptonic decays b — > clvi are excellent subjects in exploring the CKM matrix element 
V c b, we can use both the exclusive and inclusive b — > c transitions to study the CKM matrix element 
V c b- The semileptonic and nonleptonic i? c -decays have been studied extensively [8], in those studies, 
we often encounter the B c — > P, V form- factors, which are highly nonperturbative quantities and 
should be calculated by some nonperturbative theoretical approaches. In this article, we calculate 
the B* —> rj c form-factors with the three-point QCD sum rules, then take those form-factors as 
the input parameters to study the semileptonic decays B* — > rj c £Di. The QCD sum rules is a 
powerful nonperturbative theoretical tool in studying the ground state hadrons, and has given a 
lot of successful descriptions of the hadron properties [10] . There have been several works on 
the semileptonic B c -decays with the three-point QCD sum rules [HI [12l H3J [14] , while there does 
not exist work on the semileptonic B* -decays. 

The article is arranged as follows: we study the B* — > r\ c form-factors using the three-point 
QCD sum rules in Sect. 2; in Sect. 3, we present the numerical results and discussions; and Sect. 4 
is reserved for our conclusions. 

2 The B* — > rj c form- factors with QCD sum rules 

We study the B* — > rj c form- factors with the three-point correlation function Il^^^i,^), 

rV(p!,p 2 ) = i 2 [ d 4 xd 4 ye^- x e- l ^y(0\T{J 5 (x)j^O)J u (y)}\0) , (1) 
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where 

J$(x) = c(x)i~f 5 c(x) , 

^(0) = 5(0)7^(1-75)6(0), 

My) = b(yh v c(y), (2) 

the pseudoscalar current Js(x) and vector current Jy{y) interpolate the pseudoscalar meson r/ c and 
vector meson B*, respectively, the jp(0) is the transition chiral current sandwiched between the 
B* and rj c mesons. 

We can insert a complete set of intermediate hadronic states with the same quantum numbers 
as the current operators Js(ir) and J v {y) into the correlation function II^„(pi,p2) to obtain the 
hadronic representation [9] . After isolating the ground state contributions from the heavy mesons 
B* and rj c , we get the following result, 



|J 5 (0)|>? e (P2))(>? e (P2)|^(0)|g c *(Pi))(B c *(pi)|J,(0)|0) , 
(m|. -p 2 )(m 2 c -pi) 

U m tfB*m B: / . U ,2 U . A+(g 2 ) + A_(g 2 ) 

J -Wnv\rriB* +mr, c )Ai{q ) + tpi tl p2u- 



2m c (m B , - pf)(m^ - p|) ^ m B » + m Vc 

A + (q 2 ) - A_(q 2 ) a JVV) 1 

rriB*+m Va tub* + Tn Vc J 

where we have used the following definitions for the B* —¥ r\ c form-factors and decay constants of 
the B* and rj c mesons, 

A, (a 2 ) 

{VM\3^)\B*M) = ie^mB* +m nc )A 1 {q 2 ) + i{p 1 +p 2 )^e-q- 



TUB* + "V 

+iq,e ■ q A ~ { P + e^pfpl 2V <f> , (4) 
tub* + m Vc m B * + m rlc 

(0\Jl(0)\B* c ( Pl )) = f Bi m B .e^ 

f 777? 

(0\M0)\Vc(P2)) = (5) 



J2e*^(\,p)e„{\,p) = -g liV + -^-. (6) 



q^ = (pi — p 2 )^n the is the polarization vector of the B* meson and satisfies the relation 

PuPv 

X 

In this article, wg choose the tensor structures PifiP2is ? Vi^Viy 

and e^ a /3PiP2 to study the 

weak form- factors. 

Here we will take a short digression to discuss the relations among the form-factors based on 
the heavy quark symmetry [15] ■ The B* — ¥ r\ c form-factors can be re- written as 

A (a 2 ) 

(»7c(p2)|j/.(0)|flS(pi)> = iSn{mB*+m Vc )A 1 {q 2 )+i(p 1 + P 2)»e-q 2 ' 



rriB* + m Vc 

[A 3 (q 2 ) - A (q 2 )] + e^e-plp? 2V{ f ] 
q tub* + m nc 



where 



A + (q 2 ) = A 2 {q 2 ), A 3 (0) = Ao(0), 

A_(q 2 ) = -2m BS (mB S +mr,J Mq2) ~ Mq2) ■ (8) 
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Figure 1: The typical diagrams we calculate in the operator product expansion, we take into 
account the tree-level perturbative term and gluon condensates. 



In the heavy quark limit, the B* — > r\ c form-factors can be parameterized by the universal Isgur- 
wise function 

(vc(v')\j^o)\b* c (v)) = j^^Ail-v^'KH + ^/MM, (9) 

where the w M and v'^ are four-velocities, and to — v ■ v'. Then we obtain the following relations, 



V(q 2 ) = A 2 (q 2 ) = A (q 2 ) = A 1 (q 2 ) 



[m B , + m v J 2 



aco) . (10) 



The vector state \B*{v)) relates with the pseudoscalar state \B c {v)) through \B*(v)) = 2Sf\B c (v)) 7 
where the is the heavy quark spin operator. We can also express the B c — > r\ c form-factors in 
terms of the Isgur-wise function 



( Vc (v')\j,(0)\B c (v)) = £(u>)(v + v\ 



(11) 



On the other hand, the B c — > rj c form-factors are usually parameterized by the two form-factors 
F^q 2 ) and ^ (<Z 2 ), 



(vMlMVlBM) = F^q 2 ) 



(Pi +P2)^ - 



ml 



m: 



F (q 2 ) mB < 2 '^ q». (12) 



mz 



The form-factors Fi(q 2 ) and F (q 2 ) relate with the Isgur-wise function £(w) through, 



F^q 2 ) = F (q 2 



1 



(m Bt: + m nc ) 2 



mB c + m 



(13) 



Finally we obtain the following relations among the B* — > rj c and B c — > r/ c form-factors in the 
heavy quark limit, 



V{q 2 ) = A 2 (q 2 ) = A (q 2 ) = F^q 2 ) , A^q 2 ) = F {q 2 ) . 



(14) 



In the following, we briefly outline the operator product expansion for the correlation func- 
tion H fJiU (pi,p 2 ) in perturbative QCD. We contract the quark fields in the correlation function 
Yl liv {p\,P2) with Wick theorem firstly, replace the b and c quark propagators with the correspond- 
ing full propagators, then calculate all the Feynman diagrams with the Cutkosky's rule. We take 
into account the tree-level contributions and the gluon condensates in the operator product expan- 
sion, the typical Feynman diagrams are shown in Fig.l. 

For example, the lowest order tree-level diagram can be calculated as 



rV(pi,p 2 ) 



-I 



?>{ 75 [Qt+fa) + m c ] 7 „(1 - 75 ) [Qt+fc) + m b ] lu $ + m c }} 



[(k + p 2 ) 2 - ml] [(k + Pl y - m 2 } [fc 2 - m 2 ] 



ds\ds 2 



P(si,s 2 ,q 2 



(51 -P?)(*2 -Pi) 



2\ ' 



(15) 
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where 



p(si,s 2 ,q 2 



3/ 



(2tt) 3 



d 4 fc<5 [(k+ P2 ) 2 - ml] 5 [(k + Pl f - ml] 6 [k 2 - m\ 



Tr{ lb [Qt+jh) + m c ] 7^(1 - 75 ) [Q6+ fa) + m h ] + m c }} , 



(16) 



we carry out the integral over the variable k to obtain the spectral density at the quark level. The 
Feynman diagrams for the gluon condensates like the typical ones shown in Fig.l are calculated 
analogously. 

The calculations are straightforward and tedious, in the following, we present the basic formulae 
used in this article, 



d 4 kS 3 



d i kk ll S 3 



2VA(si,S2,<Z 2 ) 

7T 

2y/\( Sl ,s 2 ,q*) 



(17) 
(18) 



{aiVivViv + b 2 y 2il y 2lJ + c 2 (p lf 

iP2u H~ Ply 



where 



5 3 



5[k 2 - m 2 ]S[(k + Pl ) 2 - m 2 ]8[(k + p 2 f - m 2 2 ] 
82(81 + s 2 -q 2 ) - 2s 2 si 



bi = - 



X(si,s 2 ,q 2 ) 
Si(si +s 2 - q 2 ) - 2sis 2 
A(si,s 2 ,q 2 ) 



a 2 



b 2 = 



C2 



d 2 = 



sf+2s 2 m 2 s{sl + s 2 sf - sis 2 (si + s 2 - q 2 ) 

~" 6S2 V, , ^2 

A[Sl,S2,q ) 

S1S2 + s 2«i - SiS 2 (si + S 2 - q 2 ) 



A(si,s 2 ,g 2 ) 
3^ + 2sim 2 
\(si,s 2 ,q 2 ) 



6si 



A(si,s 2 ,(7 2 ) 2 



2sis 2 (si + s 2 - q ) - 3(sis| + s 2 sf 



T2sis 2 



si + s 2 - q" y X(si,s 2 ,q 2 ) 
Sisl + s 2 sl - sis 2 (si + s 2 — q 2 
A(si,S2,g 2 ) 2 
m 2 sxji + s 2 s1 ~ fifgCfi + s 2 - g 2 ) 
2 2A( Sl ,s 2 ,(z 2 ) 



4sis 2 



A(si,s 2 ,g 2 ) 



(20) 



^ = s 4 +m 2 — m 2 , i = l,2, and A(a, 6, c) = a 2 + 6 2 +c 2 — lab— 2bc— 2ca. The formulae in Eqs. (17-18) 
are consistent with that obtained in Refs. [161 [T7] . while the formula in Eq.(19) is slightly different 
from that of the Ref.pJ]. 

Once the analytical expressions of the correlation function at the quark level are obtained, 
then we can take quark-hadron duality below the threshold and s 2 in the channels B* and rj c 
respectively, take double Borel transform with respect to the variables P 2 = —p\ and P 2 — —p\ 
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respectively, finally obtain four QCD sum rules for the weak form-factors, 



V(q 2 ) 



2m r 



A+(q 2 ) = 



fnc m %jB*rnB*{rn B * c + m n J 

^ 3m c (si — q 2 ) + 2m b s 2 . a s GG 
12ir\(s 1 ,s 2 ,q 2 )i 7r 

2m c (rriB* + rn rlc 
fvc m njB* c m B * 



dsids2 



dsids 2 



exp ■ 
3C 



3C [m c {s l + s 2 - q 2 ) + s 2 (m b - m c 
m%, - si m 2 - s 2 ) 



M 2 



M 2 



m c + (m b - m c ) 



47rVA(si,s 2 ,g 2 ) 
si - s 2 ( Sl + q 2 + 2m 2 c - 2m 2 ) 



3siTO c 



2ir\(si,s 2 ,q 2 )?(si + s 2 - q 2 ) 

mc 

27rA(si,s 2 ,9 2 )^(si + s 2 - q 2 ) 



A(si,s 2 ,g 2 ) 



1 + 



4sis 2 



A(si,s 2 ,'? 2 ) 
Si + s 2 + 2q 2 + 



I2sis 2 q 2 



A(si,s 2 ,g 2 ) 



t a s GG ^ 



+ - 



m b 



3n\(si,s 2 ,q 2 ) 2 



1 + 



6sis 2 



A(si,s 2 ,g 2 ) 



37rA(si,s 2 ,<? 2 )2 (si + s 2 - g 2 ) 



2 , 6sis 2 (si + s 2 - q 2 ) 
Sl -2s 2 -q + T 



( ) h x P \ — 



s\ m^ c — s 2 



M 2 



M 2 



2m c (m B ' + m„J f f 3Cm c \2sis 2 - (si + m 2 - m 2 )(.S! + s 2 - q 2 )] 

— 5—7 / dsids 2 < 3 - 

/i. m t/ fl «' mfl : 1 27r 2 A(si,s 2 ,(7 2 )t 



3s 2 m c 



,a s GG 



7rA(si,s 2 ,g 2 )2 

_ TOfe 

37rA(si,s 2 ,(7 2 )i 
m c (m s . + m^J 

fvc m l c fB, 



> + 



1 + 



27rA(si,s 2 ,(7 2 )2 
6sis 2 



6sig 2 



A(si,s 2 ,g 2 ) 



A(si,s 2 ,g 2 



,« s GG.l I m 

) f ex P s — 



- m„ ) /" i 
/ dsids 2 < 



3C 



m c + (nib - m c ) 



4ir 2 v /\(s 1 ,s 2 ,q 2 ) 

4 ~s 2 (si + q 2 + 2m 2 c - 2ml) 
X(si,s 2 ,q 2 ) 



M 2 



m. 



s * , m vc- S2 
M 2 



si mi 



exp • 



S2 



M 2 



M 2 



where 



A+(q 2 ) 

A-{q 2 ) 

ds\ds 2 



C = 



A + (q 2 ) + A-(q 2 ) , 
A + (q 2 )-A_(q 2 ), 

dsi / ds 2 

(m b +m c ) 2 J Am 2 



>\2s 1 S2-(s 1 +S2-q 2 )(s 1 +ml-ml)\<y?\(s 1 ,s 2 ,q 2 )\(s 1 ,m2,ml) 



\ 



47raf 


1 


3w 


L 1 cx p( 3, J J 



Figure 2: The ladder Feynman diagrams of the Coulomb-like interactions. 



For the heavy quarkonium state B*, the relative velocity of quark movement is small, we should 

c 

account for the Coulomb-like — corrections. After taking into account all the Coulomb-like con- 
tributions shown in Fig. 2, we obtain the coefficient C to dress the quark-meson vertex, where the 
= 0.45 for the be (or 6c) systems [T3] . 

3 Numerical results and discussions 

The input parameters are taken as m Vc = 2.981 GeV [5J, m\, — 4.7 GeV, m c = 1.3 GeV, ( "* ) = 
(0.33 GeV) 4 [TBI, s i = (45± l)GeV 2 , f B . = 0.58 GeV, m B , = 6.34 GeV [19], s° 2 = (15 ± 1) GeV 2 
[9], and f Vc = 0.35 GeV [20]. In Fig. 3, we plot the weak form-factors at q 2 = with variations of 
the Borel parameters M 2 and Mf respectively. From the figure, we can see that the form-factors 
decrease monotonously with the increase of the Borel parameters at the region M\ < 2.75 GeV 2 
and Mf < 0.75 GeV 2 , and no stable QCD sum rules can be obtained. In this article, we take the 
Borel parameters as M 2 = (3.5-5.5) GeV 2 and Mf = (0.9-1.7) GeV 2 , the values are rather stable 
with variations of the Borel parameters. The contributions from high resonances and continuum 

states are greatly suppressed, exp(— -jfa) < e~ 8 and exp(— 4^) < e~ 8 2 . If we choose much larger 
Borel parameters, the numerical values of the weak form-factors changes slightly, see Fig. 3, the 
predictions still survive. The numerical values of the weak form-factors at zero momentum transfer 
are 



From Eq.(21), we can also obtain the numerical values of the weak form-factors at the squared 
momentum q 2 , then fit them to an exponential form, 



where the f(q 2 ) denote the weak form-factors Ai(0), A + (0), ^4_(0) and V(0), the C\ are c 2 are 
fitted parameters. The numerical values of the fitted parameters c\ and c 2 are presented in Table 
1. 

The calculations based on the three-point QCD sum rules indicate that the B c — > r\ c form-factor 
Fi(0) is 0.20 ± 0.02 from Ref.[TT], 0.55 ± 0.10 from Ref.JJJ], 0.66 from Ref.[13], the differences are 
rather large, as very different input parameters are taken in those studies. If the heavy quark spin 
symmetry works well, the present work indicates that Ai(0) = ^(0) = i*i(0) = 0.25, which is 
compatible with the value 0.20 from Ref.Qj]. While in the present work Ai(0) = A 2 (0) < V(0) 



MO) 

A + (0) 
A_(0) 
V(0) 



0.25 ±0.01, 
0.25 ±0.03, 
0.52 ±0.02, 
0.39 ±0.02. 



(23) 





(24) 
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Figure 3: The weak from- factors with variations of the Borel parameters M 2 and M|, where 
M 2 = 1.3 GeV 2 in (I) and M\ = 4.5 GeV 2 in (II). 

with the difference V ^y^^ — 36%, the heavy quark spin symmetry works not well enough, as 
the c quark mass is not large enough. 

The semileptonic decays B* — > r\ c ivi can be described by the effective Hamiltonian "H ff , 

H cS = ^|K b c 7Q (l- 75 )6^7 Q (l-75K, (25) 

where the V c b is the CKM matrix element and the Gf is the Fermi constant. We take into account 
the effective Hamiltonian H c g and the weak form-factors Ai(q 2 ), A + (q 2 ), A_(q 2 ) and V(q 2 ) to 
obtain the squared amplitude |T| 2 , 

\T\ 2 = 4G 2 F V 2 rve + iev«-l-vg^)( Vc (p)\j a m (26) 

where the P, p, I and v are the four-momenta of the B*, r] Cl t and P^, respectively. Finally we 
obtain the differential decay widths, 

dT = E^^ d *( P ^,p)d<%^^), (27) 

where the d$>(P —> q,p) and d§(q — > I, v) are the two-body phase factors defined analogously, for 
example, 

d$(P^q,p) = (2nfS\P-q-p) /"J . (28) 

We take the revelent parameters as G F = 1.166364 x 10" 5 GeV~ 2 , V cb = 40.6 x 10~ 3 , m e = 
0.510998928 MeV, to m = 105.6583715 Me V, m T = 1776.82 MeV from the Particle Data Group [6], 
then obtain the differential decay widths and decay widths, 

T{B* r\ c ev e ) = 5.41 x 10~ 15 GeV , 
T{B* rj c (J,y„) = 5.39 x 10" 15 GeV , 

T{B* -> t) c tv t ) = 2.13 x 10" 15 GeV, (29) 

the numerical values of the differential decay widths dT/dq 2 are shown in Fig. 4. The present 
predictions can be confronted with the experimental data at the LHCb in the future. 
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Ai(0) [ci/02] 


^+(0) [a/ci] 


A_(0) [ci/c 2 ] 


V(0) [ci/t?a] 


0.25 [0.109/0.00062] 


0.25 [0.125/0.00085] 


0.52 [0.133/0.00096] 


0.39 [0.130/0.00093] 



Table 1: The parameters for the weak form- factors. 




Figure 4: The differential decay widths with variations of the squared momentum q 2 , the A, B 
and C denote dT(B* — > ri c eD e )/dq 2 , dT(B* — > r\ c \iv ^) j dq 2 and dT(B* — > t] c tD t ) / dq 2 , respectively. 

4 Conclusion 

In this article, we study the B* — > i] c form-factors with the three-point QCD sum rules, then take 
those weak form-factors as the basic input parameters to calculate the semileptonic decay widths 
and differential decay widths. The vector B* mesons and the semileptonic decays B* — > r/ c £&(, 
have not been observed yet, the present predictions can be confronted with the experimental data 
in the future at the LHCb. 
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